On-line trading as a renewal process: Waiting time and inspection paradox 
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We briefly review our recent studies on stochastic processes modelling internet on-line trading. We 
present a way to evaluate the average waiting time between the observation of the price in financial 
markets and the next price change, especially in an on-line foreign exchange trading service for 
individual customers via the internet. The basic method of our approach depends on the so-called 
renewal- reward theorem. Assuming that the stochastic process modelling the price change is a 
renewal process, we use the theorem to calculate the average waiting time of the process. The 
so-called "inspection paradox" is discussed, which, in general, means that the average durations is 
shorter than the average waiting time. 
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I. INTRODUCTION 

Financial data have attracted a lot of attention from 
physicists as informative material to investigate the 
macroscopic behavior of markets [TH3]. Some of these 
studies are restricted to the stochastic variables of the 
price changes (returns) and most of them concern a key- 
word: Fat tails of the distributions [1 . However, also the 
distribution of time intervals can deliver useful informa- 
tion on the markets and it is worth while to investigate 
these properties extensively [3H9] and if possible, to apply 
the gained knowledge to financial engineering. 

Fluctuations in time intervals between events are not 
only peculiar to financial markets, but also very common 
in science. For instance, the spike train of a single neu- 
ron is characterized by a time series in which the time 
difference between consecutive spikes is not constant but 
fluctuates. This stochastic process, specified by the so- 
called Inter-Spike Intervals (ISI), is one of such examples 
[lOj [11]. The average of the ISI is about a few milli- 
seconds and the distribution of the durations (intervals) 
is well-described by the Gamma distribution Qj]. 

On the other hand, in financial markets, for instance, 
the time intervals of two consecutive transactions of 
BUND futures (BUND is the German word for bond) 
and BTP futures (BTP is the middle and long term Ital- 
ian Government bonds with fixed interest rates) traded 
at LIFFE (LIFFE stands for London International Fi- 
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events. 



Three typical examples with fluctuation between 





ISI 


BUND future 


Sony bank rate 


Average duration 


~ 3 [ms] 


~ 10 [s] 


- 20 [min] 


PDF 


Gamma 


Mittag-Leffler 


Weibull 



* Electronic address: j 'inoue@complex.eng.hokudai.ac.jp 
^Electronic address: Naoya.Sazuka@jp.sony.com 
^Electronic address: enrico.scalas@mfn.unipmn.it 



nancial Futures and Options Exchange) are ~ 10 seconds 
and are well-fitted by the so-called Mittag-Leffler distri- 
bution [5H7]. The Mittag-Leffler distribution behaves as 
a stretched exponential distribution in the short interval 
regime, whereas for the long interval regime, the function 
has a power-law tail. Therefore, the behavior of the dis- 
tribution described by the Mittag-Leffler function varies 
from stretched exponential to power-law at some inter- 
mediate critical interval [T2] . 

The Sony bank USD/JPY exchange rate [13], (i.e. the 
rate for individual customers of the Sony bank in their 
on-line foreign exchange trading service) is a good exam- 
ple to be checked against the Mittag-Leffler distribution. 
Actually, our results implied that the Mittag-Leffler dis- 
tribution does not well fit the Sony bank rate [15 j. The 
Sony bank rate has ~ 20 minutes [16] as the average time 
interval which is much longer than other market time in- 
tervals such as those for the BUND future. This is due 
to the fact that the the Sony back rate can be regarded 
as a so-called first-passage process [T7H22] for raw mar- 
ket data. In Table [T] we list the average time intervals 
and the probability distribution function (PDF) that de- 
scribes the data with fluctuation between the events for 
typical three examples, namely, the ISI, the BUND future 
and the Sony bank rate. From this table, an important 
question arises. Namely, how long do the customers of 
the Sony bank should wait between observing the price 



and the next price change? This type of question never 
occurs in the case of the ISI or of the BUND future, be- 
cause the average time intervals are too short to evaluate 
such informative measure. 

For the customers, an important (relevant) quantity is 
the waiting time rather than the time interval between 
consecutive rate changes. Here, the waiting time is de- 
fined as the time the customer has to wait between the 
instant in which they enter the market in the World-Wide 
Web and the next price change [13 J. If the sequence 
of time intervals is non-exponential and the customers 
observe the rate at random on the time axis, the dis- 
tribution of the waiting time no longer coincides to the 
distribution of the time intervals. 

In this review article, we present a useful way to evalu- 
ate moments of arbitrary order for the waiting time and 
for arbitrary duration distribution of price changes as 
well as observation time distribution, by directly deriv- 
ing the waiting time distribution. 

This paper is organized as follows. In the next section 
pl| we introduce the Sony bank rate |13] which is gen- 
erated from the high-frequency foreign exchange market 
rate via the rate window with width 2e yen (e = 0.1 yen 



Rate 



for the Sony bank). Then in the subsequent section III 



we explain our method to derive waiting-time distribu- 
tions. We show that our treatment reproduces the result 
obtained by the renewal-reward theorem. We also eval- 
uate the deviation around the average waiting time for 
the Weibull first-passage time distribution and uniform 
observation time distribution. We find that the resultant 
standard deviation is the same order as the average wait- 
ing time. We test our analysis for several observation- 
time distributions and calculate higher-order moments. 
The so-called "inspection paradox" [14] is mentioned in 
section IV, which means in general that the average of 
durations is shorter than the average waiting time. The 
last section \V\ contains concluding remarks. 



II. THE SONY BANK RATE: AN EXAMPLE OF 
FIRST-PASSAGE PROCESS 

The Sony bank rate we deal with in this paper is the 
rate for individual customers of the Sony bank [13] in 
their on-line foreign exchange trading service via the in- 
ternet. If the USD/JPY market rate changes by greater 
or equal to 0.1 yen, the Sony bank USD/JPY exchange 
rate is updated to the market rate. In this sense, the 
Sony bank rate can be regarded as a first passage pro- 
cesses [T7H22] . In Fig. fl] we show the mechanism of gen- 
erating the Sony bank rate from the market rate (this 
process is sometimes referred to as a first exit process 
[31]). As shown in the figure, the time difference between 
two consecutive points in the Sony bank rate becomes 
longer than the time intervals of the market rates. In 
Table [IlJ we show several data concerning the Sony bank 
USD/JPY rate vs. tick-by-tick data by Bloomberg for 
USD/JPY rate. It is a non-trivial problem to ask what 




Sony bank rate 
6 = 0.1 

e = 0.1 



FIG. 1: An illustration of generating the filtered rate by the rate 
window with width 2e from the market rate. If the market rate 
changes by a quantity greater or equal to 0.1 yen, the Sony bank 
USD/JPY exchange rate is updated to the market rate. 



TABLE II: The Sony bank USD/JPY exchange rate vs. tick- 
by-tick data for USD/JPY exchange rate. 





Sony bank rate 


tick-by-tick data 


# of data a day 


-70 


- 10,000 


The smallest price change 


0.1 yen 


0.01 yen 


Average duration 


~ 20 minutes 


~ 7 seconds 



kind of distribution is suitable to explain the distribution 
of the first-passage time. For this problem, we attempted 
to check several statistics from both the analytical and 
the empirical points of view under the assumption that 
the first-passage time might obey a Weibull distribution 
[23H25] . We found that the data are well fitted by a 
Weibull distribution. This fact means that the difference 
between successive rate changes in Sony bank fluctuates 
and has some memory. 



III. DERIVATION OF THE WAITING TIME 
DISTRIBUTION 

In this section, we derive the distribution of the wait- 
ing time for the customers. Our approach enables us to 
evaluate not only the first moment of the waiting time 
but also moments of any order. 



A. The probability distribution of the waiting time 

We first derive the probability distribution function of 
the waiting time s. Let us suppose that the difference 
between two consecutive points of the Sony bank rate 
change, namely, the first-passage time r follows the dis- 
tribution with probability density function Pw(t). Then, 
the customers observe the rate at time t (0 < t < r) that 
should be measured from the point at which the rate pre- 
viously changed. In Fig. [2j we show the relation among 
these variables r, t and s in the time axis. The waiting 
time for the customers is naturally defined as s = r — t. 



First-passage time 




Time 



FIG. 2: The relation these points r,t and s in time axis. The 
first-passage time r is given by r = ti+i —t{. The observation time 
is measured from the point t{. 



Now, notice that the distribution Q(s) can be written 
in terms of the first-passage time distribution (with den- 
sity Pw(t)) and the observation time distribution (with 
density Po(t)) of the customers as a convolution 



Q(s) oc 
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dtQ(s\r,t)P (t)P w (r). (1) 



In this equation, the conditional probability density 
Q(s|t, t) that the waiting time takes the value s provided 
that the observation time and the first-passage time were 
given as t and r, respectively, is given by 



Q(s\r,t) = 6(s -r + t) 



(2) 



where <$(•) is Dirac's delta function. Taking into account 
the normalization constant of fi(s), we have 



Q(a) 
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(3) 



where, again, t denotes the observation time for the 
customers. The result of the renewal-reward theorem 
: w = (s) = E(r 2 )/2E(r) (see for example [27]) is re- 
covered by inserting a uniformly distributed observation 
time distribution Po(t) = 1 into the above expression. 
Indeed, we have 
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where we defined the n-th moment of the first-passage 
time E(r n ) by 



E(r 
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More generally, we may choose a non-uniform Po(t). For 
this general form of the observation time distribution, the 



probability distribution of the waiting time s is given as 
follows. 



O(s) 
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where we defined 5 n by 
6. 
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By using the same method as in the derivation of the 
distribution fi(s), we easily obtain the first two moments 
of the waiting time distribution as 
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and the study of the standard deviation leads to 
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where we defined 
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Thus, this probability distribution Q(s) enables us to 
evaluate moments of any order for the waiting time. We 
consider the case of Po(t) = 1 as an example. This case 
corresponds to the result obtained by the renewal-reward 
theorem [25] . We find that S n = holds for arbitrary in- 
teger n. Thus, the waiting time distribution Q(s) leads 
to 



Q(s) = 



CPw(r) 
E(t) • 



(12) 



Then, the average waiting time and the deviation around 
the value lead to 
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For a Weibull distribution having the parameters tti, a, 
the above results can be re-written as 
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where we defined the Gamma function as 
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It is important for us to notice that for an exponential 
distribution m — 1, we have w = a = a by taking into 
account the fact that T(n) — (n — 1)!. Moreover, the 
average waiting time w is identical to the average time 
interval E(r) since w = E(r 2 )/2E(r) = E(r) holds if 
and only if m — 1 (The rate change follows a Poisson 
arrival process). These results are obtained by using a 
different method in our previous studies [25]. In Fig. |3| 
we plot the distribution Q(s) for m = 1, 2 and m = 0.585. 
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FIG. 3: The distribution of waiting time for a Weibull distributing 
Q(s) with a = 1 and m = 0.59, 1 and 2. 



B. Comparison with empirical data analysis 

It is time to compare the analytical result with that 
of the empirical data analysis. For uniform observation 
time distribution Po(t) = 1, we obtained a = 60.23 min- 
utes. On the other hand, from empirical data analysis, 
we evaluate the quantity (13) by sampling the moment 
as E(r n 



^N 



(1/7V) _C i=1 t™ directly from Sony bank rate 
data p_3] and find a = 74.35 minutes. There exists a fi- 
nite gap between the theoretical prediction and the result 
by the empirical data analysis, however, both results are 
of the same order of magnitude. The gap might become 
small if we take into account the power-law tail of the 
first-passage time distribution. In fact, we showed that 
for the average waiting time, the power-law tail makes 
the gap between the theoretical prediction and empirical 
observation smaller [T5] . 

IV. INSPECTION PARADOX 

Here we encounter the situation which is known as 
"inspection paradox". For the Weibull distribution, the 



paradox occurs for m < m c = 1. Namely, for this regime, 
we have (s) > (r) (see Fig. [4|. In general, it means that 
the average of durations (first-passage times) is shorter 
than the average waiting time. This fact is quite counter- 



<T> 



FIG. 4: Average duration (r) and average waiting time (s) as a 
function of m for a Weibull duration distribution with a = 1. The 
inspection paradox occurs for m < m c = 1. 



intuitive because the customer checks the rate at a time 
between arbitrary consecutive rate changes. This fact is 
intuitively understood as follows. When the parameter 
m is smaller than m c , the bias of the duration is larger 
than that of the exponential distribution. As a result, 
the chance for customers to check the rate within large 
intervals between consecutive price changes is more fre- 
quent than the chance they check the rate within shorter 
intervals. Then, the average waiting time can become 
longer than the average duration. 



CONCLUDING REMARKS 



As we showed in this paper, our queueing theoretical 
approach might be useful to build artificial markets such 
as the on-line trading service so as to have a suitable 
waiting time for the individual customers by controlling 
the width of the rate window. Moreover, the theoretical 
framework we provided here could predict the average 
waiting time including the deviation from empirical re- 
sults. 

We hope that this review article might help researchers 
or financial engineers when they attempt to build a suit- 
able on-line system for their customers. 
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